Theory of Weiss oscillations in the presence of small-angle impurity scattering. 
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We calculate the magnetoresistivity of a two-dimensional electron gas in the presence of a peri- 
odic potential within classical transport theory, using realistic models of impurity scattering. The 
magnetooscillations induced by geometric resonance of the cyclotron orbits in the periodic grating, 
known as Weiss oscillations, are shown to be affected strongly by the small-angle scattering processes 
dominant in conventional semiconductor heterostructures. Our results are in full agreement with 
experimental findings. 

PACS numbers: 73.20.Dx, 73.40. Kp, 73.50.Bk, 73.50.Jt 



I. INTRODUCTION 

In 1989, Weiss, von Klitzing, Ploog, and Weimann 
discovered the striking effect of a periodic potential vary- 
ing in one direction on the magnetoresistivity of a two- 
dimensional electron gas (2DEG). They observed that 
even a weak potential modulation (grating) with a wave 
vector q || e x induces very strong oscillations of the mag- 
netoresistivity p xx (B), while showing almost no effect on 
Pyy(B) and p xy (B). The oscillation period corresponds 
to a change of the ratio 2R c /a by unity, where R c is 
the cyclotron radius, and a = 2ir/q the grating wave 
length. These findings were explained by Gerhardts et 
al g] and Winkler et al || in terms of the quantum- 
mechanical band structure induced by the modulation 
(see also Q for explanation of the weak oscillations in- 
duced by the grating in p vv (B) and p xy {B) and ||||| for 
discussion of the temperature dependence of the mag- 
netoresistivity). Beenakker J(| presented a semiclassical 
theory of the effect. He showed that the geometric res- 
onance of the cyclotron motion in the grating induces 
an extra contribution to the drift velocity of the guiding 
center, with a r.m.s. amplitude oscillating with period 
2R c /a. These arguments were supported by an exact so- 
lution of the Boltzmann equation, where isotropic scat- 
tering with a relaxation time r was assumed. However, 
while accounting for the above experimental features, the 
theoretical results J||j|,|| disagree strongly with the ex- 
periment as far as the damping of the oscillations with 
decreasing magnetic field is concerned. Specifically, in 
the regime lj c t 3> 1 (which is valid for the experimental 
value of the transport time r in the whole range of the 
magnetic fields where the oscillations are observed) the 
theory predicts a linear decrease of the oscillation 

amplitude, Ap xx (B) cx B: 



Pxx (Lj%c 



(i) 



Here rj is the relative amplitude of the potential modu- 
lation ell(x) = — 7]Ep cos qx, I = vpr is the mean free 
path, w c = eB/m = vf/Rc the cyclotron frequency, Ef 



the Fermi energy, and Dp the Fermi velocity. In con- 
trast with this result, the experimental damping is much 
stronger, so that only few oscillations (6 — 7 in Q) are 
observed. 

The problem of the damping of the Weiss oscillations 
was addressed in a number of publications. In Ref. || 
an empirical observation was made that inclusion of the 
factor exp(— tt/uj c t s ) in Eq.(|l|) with r" 1 being the to- 
tal relaxation rate (as opposed to the momentum re- 
laxation rate t _1 ) yields a reasonable description of the 
experimentally observed damping. The authors of Ref. 
argued quite oppositely that only the transport time 
should enter the expression for the damping factor, and 
proposed an empirical form exp[— 47r 3 7 2 (i? c /a) 2 (w c r) _ M, 
with a fitting parameter 7. Monte Carlo simulations m, 
as well as numerical solution of the Boltzmann equation 
H , showed that allowing for anisotropic scattering leads 
to a stronger damping of the oscillations. However, de- 
spite considerable interest in this problem, no quantita- 
tive theory of the Weiss oscillations for the experimen- 
tally relevant situation, when small angle scattering is 
dominant, has been available so far. The purpose of the 
present paper is to fill in this gap. 



II. BOLTZMANN EQUATION FORMALISM 

We begin by reminding the reader briefly about the 
Boltzmann equation approach to the problem One 
starts from the kinetic equation for the distribution func- 
tion F(x, n) of electrons on the Fermi surface, 

CF(x, n) = — ew(x)En , 


C = v{x)n— +u c —-sm(j)v'{x)—-C , (2) 
Or 0(p 0<p 

where n = (cos </>, sin <fi) is the direction of the electron 
velocity, v(x) = [2m(_Ep + eU(x))] 1 / 2 is the magnitude 
of the velocity, — e is the electron charge, and C is the 
collision integral (to be discussed below). In the absence 
of the grating, the solution is given by 



1 
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(3) 



where a(°) is the Drude conductivity tensor, 
.(o) _ g o ( 1 -S 



1 + S 2 V S 1 



(To = e 2 i/^r/2 , (4) 



S = lu c t, and v is the density of states at the Fermi 
level. Using the Ansatz F — [v(x)/vf]F(°'> + one 
may transform Eq.ffl) into the form 



CF^(x,n) = V (x) v '{ x )-^{E X - SE y ) 



(5) 



This implies that the grating induced correction to the 
conductivity tensor, a = a^ + 5a, is given by 



6a = 



(Jo 



;(G(x, <j))v(x) cos <f>) 



1 -s 

S -S 2 



l + S 2 

where G(x, <p) is the solution to the equation 
CG = -2v{x)v'{x)/v 2 F = -2eU'{x)/mvj 



(6) 



(7) 



and the angular brackets (. . .) denote the averaging over 
both the velocity direction and a period of the modula- 
tion, 



(o{x,<t>))=- r dx r^o(x, 

a Jo Jo 27T 



(8) 



Calculating now the resistivity tensor p = a 1 , one finds 
that p = + 5p, where = [<7 (0) ] _1 and the only 
non-zero component of 5 p is 5p xx given by 



Spxa 



l + S 2 (Gvcosj)) 
Co 1 + (Gv cos (f)) 



(9) 



For a weak modulation, r] -C 1, the corrections dp and 
5a are proportional to rj 2 . Therefore, to calculate the av- 
erage (Gv cos 4>) in the leading order, one would have to 
solve the Boltzmann equation (Q) up to the second order 
in r\. However, it was observed by Beenakker || that the 
following identity holds: 



(G(x, 4>)v(x) cos (p) 



CT 



m{l + S 2 ) 



(G(x,cp)U'(x)) . (10) 



Since U(x) cx rj, it is now sufficient to calculate G(x, </>) to 
first order in 77, in order to find 5a or 5p to the order rj 2 . 
We note that, while only presented for the isotropic scat- 
tering case in Ref . H , the identity ( |l0| ) holds for a general 
collision integral C as well, since its derivation uses only 
the fact that C{1} = 0, C{cos0} = — (1/r) cos</>, and 
C{sm<fi} = — (1/r) sin</>. 

To proceed further, one has to solve the kinetic equa- 
tion (^). This requires first specifying the collision inte- 
gral. Its general form reads: 



C{F}(<j>) = / d<t>'P& - cf>')[F(ct>') - F(4>)\ , (11) 



where P(<f> — <fr') is the scattering cross-section. The col- 
lision integral can be most conveniently characterized by 
its eigenvalues 

C{e"™*} = -A m e im * ; m = 0, ±1, ±2, . . . 

A m = X-m = / d<j)P(<p)(l - cos 7720) > . (12) 



Note that Ao = reflects the particle number conser- 
vation; Ai = 1/r is the transport relaxation rate, while 
Aoc = 1/t s is the total relaxation rate. Expanding the 
solution of Eq.(^|) in the angular harmonics, G(x,<f>) = 
Sm=-oo Sme'^e'^ + c.c, we get to the leading order in 
r\ the following coupled system of equations PQ |: 



0>n9n — 9n+l + flVi-1 + C n J 

2l 77 

a n = (tnuj c + A n ) ; c n = 6 n0 

VFq vf 



(13) 



According to Eqs.(|), ©, the correction to the resistiv- 
ity tensor is determined by go , since 



^\ 779-Eft 

(Gv cos 0) = — — — j-Imgo. 

m(l + 5 Z ) 

In turn, go is found from the system (|0|) as 
77 1 



,9o 



7JF i?l + R-i 



(14) 



(15) 



where R\ = 51/50 is determined by the set of (homoge- 
neous) equations ( |l3| ) with n > 1, and = i?i| Wc -,- Wc . 
For a generic set of eigenvalues A„, the solution R\ can 
be presented as a continuous fraction, 



Ri = 



1 



a 1 



(16) 



The result can be expressed in closed form in the case 
when A n depends linearly on n, X n = t _1 [1 + 2p(|n| — 1)] 
with an arbitrary p, yielding |Io[| 



R\ — —J\-/J- ; R-i — J\+/J+ , 
J± = J±(i -2p)i/(S±2pi) (Q/(l±2pi/S)) ; 

Jl± = Jl±(l-2p)i/(S±2pi) 

(Q/(l ± 2pi/S)) 



(17) 



where Q = qR c - The two most important particular 
cases are p = 0, corresponding to a white-noise random 
potential (isotropic scattering), and p = 1, describing 
the scattering by a white-noise random magnetic field 
with the scattering rate P(4>) = t _1 cot 2 (0/2). In the 
next section we investigate 5p xx in these two situations, 
exploiting the existence of the exact analytical solution 
( |l7| ) . This will give us a first example of how the small 
angle scattering (dominant in the random magnetic field 
case, when r" 1 = 00) modifies the magnetic field depen- 
dence of the amplitude of the Weiss oscillations. 
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III. WHITE-NOISE RANDOM POTENTIAL 
VERSUS RANDOM MAGNETIC FIELD 

Using Eqs. Q), (|l5|), and we find for arbitrary 
values of the parameter p 



(Gwcos0) = 



rj 2 ql 



2(1 + S 2 ) 



J 



For the particular case of isotropic scattering, p = 0, the 
following identity can be used to simplify Eq. : 



J_Ji_ 



7 / - 21 
J+Jl - ~ SQ 



sinh(7r/S') 



tt/S 



This allows us to reduce Eq.(|18[) to the form 

v 2 (qi) 2 J+J- 



{Gv cos <fi) — 



4(1 + 52) sinhfr/g) 

71"/ S 



J+J- 



(19) 



(20) 



Substituting © in (|) and assuming that 77 is small 
enough so that (Gv cos ^) < 1 (this will be valid for all 
B, if rj 2 ql/4: < 1), we finally get 



pQ 4 W sinh(7r/S) 



2 Jj/s{Q)J-i/s(Q) 
22^-J 4/s (Q)J_ l/s (Q) ' 



(21) 



where po = o - ^ 1 is the Drude resistivity. Eq.(^i"|) is com- 
pletely equivalent to Beenakker's result (Eq. (5) of Rcf. 
Q), which is, however, presented in Q in the somewhat 
different form of an infinite series. The equivalence fol- 
lows from the Bessel function identity lilt] 



E 



(22) 



We turn now to the random magnetic field scattering, 
where p = 1. Analyzing Eq. ( |l8| ) in this case, we find 
the following two asymptotic regimes: 



$Pxx V 



Po 



^-S 2 cos 2 (Q - tt/4) , 4Q/5<1 



s 



(25) 



i2Q] 



4Q/S> 1 (26) 



Again Eq.(p5|) corresponds to a power-law behavior of 
the oscillation amplitude in the range of relatively high 
magnetic fields B, while Eq.p^) describes the exponen- 
tial damping of the oscillations in low fields. There is 
however a number of important differences between the 
results for the isotropic scattering [Eqs. (23), ( Ejj) ] and 
those for the random magnetic field scattering [Eqs.(p5|), 



i) the "high" field behavior ( |25|) is Spxx/po oc B 2 7 
rather than Sp X x/po oc B as in Eq.(|23[); 

ii) the value of the magnetic field, at which the expo- 
nential damping starts, is much larger in the ran- 
dom magnetic field case (oj c t ~ 2y/ql ^> 1), than 
for the isotropic scattering (lu c t ~ 7r); 

iii) the low-field damping is of a Gaussian form (with 
respect to 1/B) in the random magnetic field case, 
{dpxx/ Po)osc oc exp[— 4ql/(cu c T) 2 ], rather than of a 
simple exponential, exp(— tt/uj c t), as in Eq.d24). 



All these features are consequences of the small- angle 
scattering, dominant in the random magnetic field case. 
In the next two sections we will consider the Weiss os- 
cillations in the long-range correlated random potential, 
which also leads to small-angle scattering. Since in this 
case the exact solution can not be expressed in a simple 
analytic form, we will use different methods to find the 
shape of the Weiss oscillations. 



We will assume that ql = QS 3> 1, which is a nec- 
essary condition for the existence of oscillations. At 
Q = qR c 1 (i.e. in the oscillation region), the asymp- 
totic behavior of Sp xx /po reads according to Eq. ( pl| ) 

*M = "f-QS 2 cos 2 (Q - tt/4) , n/S « 1 (23) 

PO 27T 



5px 



po 



QS[l + 2e-^ /s sm2Q] , n/S > 1 (24) 



Thus, the amplitude of oscillations depends linearly on 
magnetic field B in the high- field regime ui c 3> tt/t, 
and gets exponentially damped in 1/B, (Sp xx / p ) osc oc 
e -7r/w c r^ m i owes t fields, uj c <C tt/t. Note that for 
a typical value of the transport time r in the GaAs het- 
erostructures, r ~ lOOps, the condition tt/lo c t ~ 1 corre- 
sponds to very low fields, B ~ 0.01T. 



IV. LONG-RANGE RANDOM POTENTIAL: 
EXPONENTIAL DAMPING OF WEISS 
OSCILLATIONS 

The random potential experienced by electrons of a 
2DEG in the GaAs heterostructures is very smooth, since 
the impurities are separated from the 2DEG layer by a 
relatively large distance (spacer) d s > kp 1 (&f being 
the Fermi wave vector). Assuming that the positions of 
ionized donors are uncorrelated, while their density is 
equal to that of the electrons, n e = k F /2ir, and treat- 
ing the screening in the random phase approximation, 
one finds for the random potential correlation function 
W(\r — r'|) = (U(r)U(r')) in momentum space 



W(q) = 



{2vf 



(27) 
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where v — m/2ir is the density of states per spin orien- 
tation. The scattering cross-section in the Born approx- 
imation is thus given by 

P{6) = vW(2k F sin 0/2) = ^%- 4fe ^ sin . (28) 

2m 

The total (t^) and the momentum (t _1 ) relaxation 
rates are defined as 



l/r s 



2 7T 



4meL 
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1/r = / d4>P{(j)){l - cose 
j o 



1 



(29) 
(30) 



so that their ratio is r/r^ = (2kpd s ) 2 ^S> 1. Both t and 
r s are experimentally measurable: the former is found 
from the dc conductivity at B = 0, while the latter 
follows from the Dingle analysis of the damping of the 
Shubnikov-de Haas oscillations. For the structures of 
the type used in the Weiss oscillation experiment Jl]] 
(n e ~ 3 x 10 11 cm -2 , d s ~ 30nm), their typical values 
are r ~ 50ps, r s ~ 3ps, so that the ratio t/t s ~ 20 is 
large, in agreement with the theory. Note, however, that 
the experimental values of r and r s are usually some- 
what different from the theoretical estimates (|29|), (30). 
This was attributed to correlations between the positions 
of the scatterers |f2],|l3). Below we will express our re- 
sults for the damping of the Weiss oscillations in terms 
of r and t s , so that for comparison with experiment one 
should simply use the experimental values of these quan- 
tities. 

As was shown in || , the physical origin of the oscilla- 
tions can be traced to the additional contribution to the 
drift velocity of the guiding center, 



1 



t'dr 



2ttB 



d(j)E(X + R c cos< 



(31) 



where X is the ^-coordinate of the guiding center and 
E(x) — —U'(x) = — (i]EFq/e) s'mqx is the electric field 
induced by the grating. Averaging v\ r over X , one finds 



2 \2-KeB J jq 



r> 2 e 



iqR c (sin 0i— sin $2) 



(32) 



For qR c ^> 1, the oscillations are determined by the in- 
terference of the contributions from the vicinities of the 
saddle points <j)\ = 0, <j>2 = n (and vice versa), leading to 



J drlosc 



cx cos(2qR c — n/2) 



(33) 



(in this section we do not keep power-law prefactors, since 
we calculate the exponential damping factor only). 

Impurity scattering leads to a change AX of the guid- 
ing center position X during the time an electron moves 



along the cyclotron orbit from (f> = to 
introduces in Eq.(p3[) a damping factor 



D 



JqAX\ 



7r. This 



(34) 



For isotropic scattering, any scattering event kicks an 
electron far away from its original cyclotron orbit, so 
that D is determined simply by the probability for an 
electron to remain unscattered during the half cyclotron 



period time n/u c , yielding D ~ e 



-7T /iU c T 



in agreement 



with Eq.(|24j). The situation is, however, different when 
the scattering is of small-angle character, so that even 
after one or several scattering events an electron remains 
close to its original cyclotron orbit and contributes to the 
oscillations. 

We first consider the random magnetic field scatter- 
ing. Since in this case the exact solution has been found 
(Eq. (^6|) , this will be an additional check for the correct- 
ness of our method of calculation in this section. 

In this case, the probability for an electron to be scat- 
tered by an angle 5<p in a small time interval St is given 

by @ 



V{84>) 



1 



2St/r 



ir (2St/r) 2 + (Scj)) 2 



(35) 



Scattering by the angle 8(f) changes the position X of the 
guiding center by R c sin 4>8<f>, so that 



AX = R c 



k 



where 4> k = <t>{tk) = u c t k \ t k = kir/Nuj c 
N — > oo. According to Eq.(B9), 



1,2, 



d{5(j>)V{8(l>y z5 <t' = e - 2zSt / T . 



Therefore, the damping factor D is equal to 

JJ _ TJ^ e «9-R<: sin 0fe<50fe^ _ J^J g-2g-Rc sin <p k &t k /r 
k k 



(36) 



(37) 



j 2qR c r . \ f AqR c 
exp < / dfflsin© > — exp < > , (38) 

/O J I U C T 



U! r T 



in precise agreement with the exact solution, Eq.(|26|). 

We turn now to the long-range random potential. In 
this case, the scattering cross-section ( p8| ) implies that 
the distribution function of the scattering angles 8<p in a 
time interval St is given by 



ruM = ( i- 5 1 ) ()((), >) + si r(<u >) 



Therefore, 



(e lz54 ') = 1 - 



St 



t s z 2 + {2k F d s ) 2 



(39) 



(40) 
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This leads to the following result for the damping factor 



,iqR c sin <p k 84> k \ 



,2 p 7T 



exp 



exp 



(qRc 

^cTs Jo (qR c sin 4>) 2 + (2kpd s ) 2 



1 - 



^/l + (T s /T)(qR c y 



(41) 



As is seen from Eq.(|4l]), the region of exponentially 
damped oscillations can be further subdivided into two 
regimes: D ~ exp{— n(qR c ) 2 /2uj c t} at (qR c ) 2 <C t/t s 
and D ~ exp{— it/uj c t s } at (qR c ) 2 3> t/t s . In the for- 
mer the leading contribution to D comes from multi- 
ply scattered carriers, while in the latter D is governed 
by those electrons, which remain unscattered within the 
time tt/u) c . 

The exponential damping found in the above is efficient 
at low enough magnetic fields, uj c t < [(tt /2)(ql) 2 ] 1 / z . At 
stronger magnetic fields, the independence of Sp xx is of 
a power-law nature. This will be considered in the next 
section using a different method. 



V. LONG-RANGE RANDOM POTENTIAL: HIGH 
MAGNETIC FIELDS 

To find the grating-induced correction to the resistiv- 
ity at relatively strong fields, we return to the Boltzmann 
equation approach of Secjl|. For the scattering cross- 
section (|2^) the eigenvalues A/ of the collision integral 
are equal to 



A; = 7r — / d(pe 
m Jo 



I 2 



4feFdsSin «V 2 (l-cos^) 

COS l(j>) 



{t s /t)1 2 + 1 



(42) 



The approximation we will make at this point will be to 
consider the limit t s /t — > 0, which reduces Eq.(f42|) to 
A; = I 2 jr. This approximation is valid for not too small 
magnetic fields, such that lu c qvp^s/r) 1 / 2 , where the 
total relaxation rate r s is irrelevant, and only the trans- 
port one, r, is important. This will allow us to find 
5p X x in the high-field regime (where the independence 
of the oscillation amplitude is of the power-law nature) , 
which was outside the regime of validity of the consid- 
eration in Sec.[j"v|, where only the exponential damping 
factor was studied. Furthermore, we will also repro- 
duce the result of Sec. IV for the intermediate regime 
(2/tt)oj c t <C (qR c ) 2 "C t/t s , where the damping is al- 
ready exponential, but r s is still irrelevant. 

Within the above approximation, the collision in- 
tegral C{G} is replaced by the differential operator 



r- 1 d 2 G/d 2 (j) 2 , and Eq.(0) with G(x, 0) = g(</>)e iqx + c.c. 
takes the form 



^ + iqR c cos < 
acp 



1 d 2 q q 

— /d? = iT iT-- ( 43 ) 

Ul r T OCT ILOr 



For lo c t 3> 1, ql 3> 1, we can solve this equation itera- 
tively considering the d 2 g/dcf> 2 term as a perturbation. 
Dropping this term in (43) and differentiating with re- 
spect to 4> we find 



d 2 g 



iqR c g sin 



dg 



■ COS( 



(44) 



Substituting this back to Eq. (ff3[) , we get a first order 
differential equation, the solution of which reads 



<?(0) 



iqrj 
2lo c 

x exp 



Q2 

-iQ(sin0-sin0') - — ( 



(45) 



Here we introduced again the compact notations Q — 
qR c , S — ui c t. Substituting now the solution d45| ) into 
Eqs. (g), (|l5|), we finally get 



Spxx _ V 2 g2 irQ 2 /2S 

~^o~ ~ Y S inh(nQ 2 /2S) JiQ2/2S 



(Q)J- 



iQ 2 /2S 



(Q) 



(46) 



Combining the results of this and the preceding section, 
we get a complete description of 5p xx in the full range of 
magnetic fields B, which will be analyzed in Sec. VI. 



VI. LONG-RANGE RANDOM POTENTIAL: 
ANALYSIS OF THE RESULTS 

According to the results of Sees. |lV|, [y], we find in gen- 
eral the following three regions of the magnetic field B 
with different behavior of the oscillation amplitude: 

i) high field domain, Q < Q x , with Qi = (2ql/ir) 1 / 3 . 
In this region irQ 2 /2S <C 1, and Eq.(E|) yields 



Sp x 



Po 



?? 2 -^-cos 2 (Q-7r/4) 



(47) 



The amplitude of the oscillations is proportional to 
B 3 in this regime, in contrast to the linear behav- 
ior, Eq.(p3|), for the isotropic random potential, and 
to the B 2 behavior, Eq.(p5|) in the case of random 
magnetic field scattering. 

ii) intermediate domain Qi <C Q <C Q2, with 
Q2 = -\/2t/t s . The oscillations get exponentially 
damped: 

^ = ?F-QS[1 + 2e-^ 2 ' 2S sin 20] . (48) 
Po 4 



5 



iii) low- field region, Q 3> Q2- The exponential damp- 
ing factor changes its form: 



Po 



= ^QS[1 + 2e-*l s > sin 2Q] , (49) 



where S s = oj c t s - Note that Eq.(|49|) differs from 
the low-field behavior for the isotropic scattering 
case, Eq.(^4j), only in that r is replaced by r s in 
the exponential damping factor. 



The two regimes (i) and (ii) are jointly described by 



Po 



-Q 2 S 



sinh 7T[i 



(50) 



with /i = Q 2 /2S, while both the regimes (ii) and (iii) are 
given by 



Spxx 

Po 



n 2 

\QS 



l + 2exp<J - — 



T s n\ -1/2' 

1 + -Q 

T 



sin 2Q 



(51) 



Finally all the three regimes can be described by a single 
formula of the form d5Q) with 



1 



1 



1 



T \ -1/2' 
T 



(52) 



Let us also note that the behavior of 8p xx at Q -C 1 
(i.e for magnetic fields stronger than those generating the 
oscillations) is the same for all the types of scattering, 



8px 



n 2 

'/ Q 2 

T 



(53) 



as is easily found from Eqs. (|l8|), (pi|), (|4q). In reality, 
in most of the Weiss oscillation experiments the magne- 
toresistivity at such strong fields is strongly affected by 
the Shubnikov-de Haas oscillations and by the quantum 
Hall effect. However, the giant quadratic magnetoresis- 
tance (^3|) was observed in the experiment |lq| , where 
the measurements were performed at somewhat higher 
temperatures. 

Now we evaluate the boundaries Qi, Qi of the re- 
gions with different behavior of the oscillation amplitude 
in a typical experimental situation. We find from the 
parameters of the sample of Ref. M = 27r/382nm, 
n e = 3.16 x 10 11 cm -2 , I = 12/jm) the transport time 
t = 52ps. Furthermore, a typical value of the total re- 
laxation rate for such structures is ]i"2|] r s ~ 3ps (as we 
will see below, with this value of r s our results perfectly 
describe the experimentally observed damping of the os- 
cillations in low fields) . This gives the following crossover 
values: Q± ~ 5.0, Q2 — 6.0. For a higher mobility sam- 
ple studied in Ref. |l6| (the Weiss oscillation data are 
presented there in the inset of Fig. 1) these values are 



somewhat larger: Qi ~ 8.5, Q2 — 11- We conclude that 
although parametrically Q2/Q1 oc (fc^/g) 1 / 3 is large at 
kp S> q, for realistic parameters Q2 and Qi are numeri- 
cally close to each other. Therefore, the regime (ii) corre- 
sponds in reality to a rather narrow intermediate domain 
between the regime (i) describing the fast (cx B 3 ) drop 
of bp xx in high fields, and the regime (iii) corresponding 
to the exp(— tt/lu c t s ) damping of the oscillatory part of 
bp xx in low fields. 

In Fig.l we compare our theoretical results with exper- 
imental data of Weiss et al . The sample parameters 
are @ q = 27r/382nm, n e = 3.16 x 10 n cm~ 2 , r = 52ps, 
the total relaxation rate is taken to be r^ 1 = (3ps) _1 . 
As is seen from the figure, with the modulation strength 
7/ = 0.065 we get a very good description of the experi- 
mentally observed magnetoresistivity [fl7| . 
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FIG. 1. Experimental data of Ref. |jj (dash-dotted line) 
compared with our theoretical results for the long-range po- 
tential scattering (full line). Parameters are: q = 27r/382nm, 
n e = 3.16 x 10 n cm- 2 , r = 52ps, r s = 3ps, rj = 0.065 



To illustrate the difference between the long-range po- 
tential scattering and the isotropic scattering, we plotted 
in Fig. 2 the corresponding results for the modulation- 
induced resistivity correction at the same value of the 
transport time r. It is seen that the results for these two 
different scattering mechanism differ drastically not only 
in low magnetic fields (where the small-angle scattering 
leads to the exponential damping of oscillations), but also 
in high fields (due to the different power-law dependence 
of 8p xx on B). In particular, already at the first max- 
imum of Sp xx , the value of the correction is smaller by 
factor of a ~ 5 for the long-range random potential scat- 
tering, as compared to the isotropic scattering 
the same value of the grating strength 77. For the second 
maximum this ratio is already a factor of ~ 10, etc. As 
a consequence, one generates a considerable error when 



6 



trying to evaluate the modulation strength by fitting the 
amplitude of one of the first maxima by the white-noise 
random potential results. For example, Beenakker ||] 
underestimated r\ by factor of 3 plf . 



8000.0 




B (Tesla) 



FIG. 2. Grating-induced correction to the 2DEG reistivity, 
dpxx/po, in units of r/ 2 /4 for the isotropic potential scatter- 
ing (r = 52ps, full line) and the long-range random potential 
scattering (r = 52ps, t s = 3ps, dash-dotted line). The sample 
parameters (q — 27r/382nm, n e = 3.16 x 10 11 cm -2 ) are the 
same as in Ref. |l| and in Fig.l. 



magnetic fields, the oscillatory part of the magnetoresis- 
tivity is found to be exponentially damped, much more 
efficiently than in the case of isotropic scattering with 
the same r. The low-field damping factor is of the form 
e -7r/u c T,^ familiar from the theory of the Shubnikov-de 
Haas oscillations. Let us stress, however, that while the 
Shubnikov-de Haas oscillations are of purely quantum na- 
ture, the Weiss oscillations represent a classical effect. 
Nevertheless, due to the interplay of the finite modula- 
tion wave vector q and the cyclotron radius R c , the resis- 
tivity correction Sp xx is determined by the whole spec- 
trum of eigenvalues A; of the collision integral (and, in 
particular, by Aoo = 1/t s ), rather than by Ai = 1/r only 
(governing conventional transport in the absence of mod- 
ulation) . Finally, our results describe nicely the shape of 
the experimentally observed magnetoresistance and the 
dependence of the oscillation amplitude on the magnetic 
field. 

We are grateful to D. E. Khmelnitskii for attracting 
our attention to the problem of Weiss oscillations and 
to Y. Levinson for useful discussions at various stages of 
this work. This work was supported by the SFB195 der 
Deutschen Forschungsgemeinschaft . 



VII. CONCLUSIONS 

In this paper, we have presented a theory of the Weiss 
oscillations of the magnetoresistivity in a 2DEG, for the 
case when the impurity scattering is predominantly of 
the small-angle nature. While this is precisely the exper- 
imental situation in currently fabricated semiconductor 
heterostructures, no quantitative theoretical description 
of this case has been available so far. Previous theo- 
ries assumed the scattering to be isotropic. We 
have shown that the small-angle scattering strongly af- 
fects the shape and the damping of the Weiss oscilla- 
tions. This was demonstrated both for the scattering 
by a white-noise random magnetic field (linearly increas- 
ing eigenvalues A; = t _1 (2|Z| — 1) of the collision in- 
tegral) and for long-range random potential scattering 
(A; ~ t -1 P /[(t s /t)1 2 + 1]). In the latter case, which was 
of primary interest for us, the shape of the magnetoresis- 
tivity is determined by both the momentum and the total 
relaxation rates (t _1 and TjT 1 , respectively), and has the 
following features. In relatively strong magnetic fields B, 
the amplitude of the modulation induced correction Sp xx 
falls of as B 3 with decreasing B, i.e. much more rapidly, 
than in a white-noise random potential, where Sp xx oc B. 
This rapid fall-off is determined by the quadratic increase 
of the eigenvalues A; oc I 2 for not too large /. In low 
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